Expanding, axisymmetric pure radiation gravitational fields with simple twist 
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New expanding, axisymmetric pure radiation solutions are found, exploiting the analogy with 
the Euler-Darboux equation for aligned colliding plane waves. 
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I. INTRODUCTION 

o 

■ There exist many papers dealing with algebraically special, expanding and twisting pure radiation solutions of the 
Pn| I Einstein equations [1-7]. The standard form of the metric is [1]: 

^' ds^ = ?^-2n(dr + WdC + WdC + Hn) (1) 
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n:^du + LdC + LdC (2) 

Here r is the coordinate along the null congruence of geodesies, u is the retarded time, the complex coordinates C, C 
span a two-dimensional surface. The metric components are determined by the r-independent real functions P, m, 
M and the complex function L: 

P = — ^ (4) 



g^; W ^ p-^Lu + idY. (5) 
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W). H = -r (In P)„ - (mr + MT^pp + — (6) 



K = 2P'^Re [d (5 In P - L„)] (7) 



' where 9 = 9^— Ldu and E is the twist. The basic functions P, L, to, M satisfy the system: 

(<9 - 3i„) (to + iM) = (8) 

p-^M = ImddddV (9) 

^ -2P^ [P-^ (to + iM)]^ + (^ddddV)^ - 2P^ (ddV)^ (10) 

where Vu = P, n is the energy density of pure radiation and the Newton constant is set to 1. Eqs. (8), (9) and (10) 
are in fact Eqs. (26.32) and (26.33) from ref. [1]. 

It has been noticed in different contexts that the condition M = (vanishing NUT parameter) simplifies the 
equations [4,5,6,8,9]. Twisting gravitational fields with M = generalize the classes of Robinson- Trautman [10] and 
Kerr-Schild fields [11] which are physically realistic, their simplest representatives being the Schwarzschild, Kerr and 
Vaidya solutions. 

In the present paper we explore this condition applying the method of Stephani [2] . We discuss axisymmetric fields 
with the simplest possible twist. In Sec. II Eqs. (8-10) are reformulated in terms of an invariant potential which 
leads to the i„ = gauge. In Sec. Ill the main equation (9) for simplest twist is shown to be equivalent to the 
Euler-Darboux equation, which is central in the theory of aligned colliding plane waves (CPW). We use the known 
solutions and techniques to find solutions for our problem. In Sec. IV a closing discussion is presented. 
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II. FIELD EQUATIONS IN THE = GAUGE 

Following [2] we introduce the invariant complex potential 4> which solves Eq. (8): 

m + iM = <t>l (11) 

-I 

When M = we can apply the gauge transformation 

u' = f{u,a) (13) 

(m + iM)' = (m + iM) (14) 
to make the mass parameter m a positive or negative constant mo so that 

(f) = [u + iq (a)] (15) 

L = iCg« (16) 

where q is real and due to the axial symmetry depends only on ct = CC The complex coordinate Q is related to the 
angular coordinates 0, on the distorted spheres {r = ro, u = uq ) according to 

C = \/2 tan (61/2) e'^ (17) 

Obviously i„ = 0. This gauge differs from the usual Kerr's gauge [12] P„ = 0, but is very suitable when the NUT 
parameter vanishes. Eqs. (9) and (10) simplify 

ddddV = ddddV (18) 

= GmoP-^Pu + 2P^ddddP - 2P^ddPddP (19) 

The second equation is in fact an inequality. When P„ ^ 0, can be made positive by the choice of mo at least for 
some region of spacetime [1,4,8]. The expressions for the metric components simplify too, e.g. the gauge invariants S 
and K read 

S = (20) 

K = P^ {dd + dd)lnP (21) 

where Q = q^-^. 

When m + iM = (Petrov types III and N) Eq. (8) is an identity but still a potential (f> may be introduced with 
the property d(j> = and the subclass of solutions satisfying Eqs. (15) and (16) (with mo — 1 ) can be studied. This 
results in setting mo = in all other equations. 

In both cases the main equation (18), which is of fourth order with respect to V, becomes a linear second order 
equation for P. Let us choose the simplest possible twist q = a, Q = 1, L = i(. Then Eqs (18) and (19) read 

{dd + dd)P = (22) 



= 6moP"^P„ - GP^Puu - 2a^P^ 



2P„„ + -P<, +4P 



(23) 



The last term in Eq. (23) is obviously negative, so necessarily the first must be positive for a type II solution and the 
second must be positive for a type III or N solution. 
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III. REDUCTION TO THE EULER-DARBOUX EQUATION 



Let us introduce the complex variable z = ^ (u + m). Then Eq. (22) becomes 

2 (z + z) P,,- + + = (24) 

When z and z are two real variables this is the Euler-Darboux equation, the main equation in the theory of aligned 
CPW [13]. We can adapt the numerous solutions for our problem. We must ensure that P is real and investigate the 
regions where > 0. 

In the original variables Eq. (24) may be written as 

Pun + Paa + -Pa = (25) 

a 

which is the analog of the equation for vacuum Gowdy cosmologies [14]. Eq. (25) has been derived from another 
viewpoint in [8] and discussed there. Solutions with separated variables behave like (m + c) In ct or e""Jo {ha) where 
a, 6, c are constants and Jq is a Bcsscl fimction. The first class includes the simplest CPW solution P = — In (2 + 2). 
A simple solution is obtained by separating the variables in Eq. (24) [13]: 

P = A{{a^z){a-z)\-^l^ (26) 
where A, a are constants. A is ignorable, while a must vanish for a real P. Hence 

P = 5-^2 (27) 

B = (7^ + (28) 

The energy density becomes 



n 
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-6mouB"^ - 12B-3 (29) 



We shall discuss only positive retarded times w > uq > 0. The energy density is regular in u and a. If mo < and 
—mo > 2/mq then r? > 0. Type III solutions, however, always have negative energy and are not realistic. The gauge 
invariants are also regular and vanish for big u 

S = B-^ (30) 

K = -2(75-2 (-3^^ 
The Weyl scalars [1,15] have the following leading terms 

*2 = mop^ (32) 

*3 = -p'P'^dl + O (33) 

*4 = pP'/„ + 0(p') (34) 
where / = P'^ddP. Plugging Eq. (27) into Eqs. (33) and (34) we obtain 

*3 = -6p2^(a + m)'S-7/2^0(p3) (35) 

*4 = QipC {o- + iuf B-^ + O (p^) (36) 

It can be shown, using the full expressions for '^^ and ^4 [15], that the Weyl scalars are regular in a and u. When 
u 00, ^"3 and \l/4 vanish, while \l/2 —mo/r^. 

There is an analog of the cosh"^ solution for CPW. Suppose that P = P{y) where 
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.= ^t^ = - (37) 

z + z a 

Then Eq. (24) yields the solution 

P = sinh"^y (38) 

The energy density is given by the expression 

^2 ^ 6mo ^ 6P3,/ ^ 2p2y2 (t/^ + 4) 

MP(l+y-2)l/2 y2(y2 + 1)3/2 ^2(^2 + 1) 

The first term is regular in y and positive when mo > 0. When y — > oo the energy density becomes negative and 
^ _y2 jjj2 ^ matter how big mo may be. 

Eq. (25) is the starting point for the procedure leading to the first Yurtsever solution [13,16]. It comprises the 
expressions 

P = B^/^Pi {x) (40) 

P = B^'^Qi {x) (41) 

where Pi and Qi are Legendre functions of the first and second kinds and x = uB~^/'^. Solutions with P; grow 
infinitely when tr — > oo, u — > oo and change sign. A promising solution is 

P = 2Qo(x) =lnii| (42) 

It's energy density is 



The range of x is < x < 1 and the first two terms arc always positive when mo > 0. The last term has a negative 
pole for cr ^ 0, ( a; ^ 1 ) of the type cr"^ ^ can't be compensated by the positive singularity of the second 

term which is ~ In^ cr. Hence < for cr ^ 0. 

In [13] a general method is presented for obtaining solutions of the Euler-Darboux equation. Angular coordinates 
arc introduced, which is possible for CPW because the variables are bounded. This method docs not have an analog 
for expanding and twisting solutions. 

At last, it should be mentioned that any linear combination of the solutions derived above is also a solution of Eq. 
(24). 



IV. DISCUSSION 



We have shown that when the NUT parameter vanishes and the gauge L„ = is used, the main equation (22) for 
axisymmetric expanding pure radiation fields with the simplest twist becomes the second order, linear, Euler-Darboux 
equation for P. This is the central equation in the theory of aligned colliding plane waves. We have found analogs of 
some of the numerous known solutions adapted to our problem and studied the region of spacetime where the energy 
density of pure radiation is positive. 

It is interesting that the interaction region of two colliding aligned plane waves is mathematically equivalent to an 
expanding distorted spherical wave with a simple twist. However, this is not the first example of such kind. It was 
shown in [17] that the field equations for non-aligned colliding electro-gravitational plane waves coincide with the Ernst 
equations [18] for stationary axisymmetric Einstein-Maxwell fields. The Bell-Szekeres interacting solution is equivalent 
to the twistless conformally flat Bertotti- Robinson solution [13] . Colliding gravitational waves are algebraically general 
due to the presence of both ^'i and 5*4. When an electromagnetic wave collides with a gravitational one sometimes the 
solution is of type II [19] without twist and falls in a class of solutions generalizing the Robinson- Trautman solutions 
(distorted spherical twistless waves) [20] . 

The structure of time-dependent expanding and twisting algebraically special pure radiation solutions obviously is 
very rich, because in the case of simplest possible twist they coincide mathematically with the large variety of aligned 
colliding plane waves. 
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